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Abstract. We report on an observation of the recently discovered accreting millisecond X-ray pulsar IGR J0029 1+5934 per- 
formed with the i^Xr^-Proportional Counter Array (PCA) and Chandra-High Energy Transmission Grating Spectrometer 
(HETGS). The RXTE data are from a two week follow-up of the source while the Chandra observation took place around 
the end of the follow-up, about 12 days after the discovery of the source, when the source flux had decreased already by a 
factor of ten. The analysis of the Chandra data allowed us to extract the most precise X-ray position of IGR J0029 1+5934, 
RA=00'' 29'" 03.08' and Dec=+59° 34' 19.2" (0.6" error), compatible with the optical and radio ones. We find that the spec- 
tra of IGR J0029 1+5934 can be described by a combination of a thermal component and a power-law. Along the outburst 
detected by PCA, the power-law photon index shows no particular trend while the thermal component (~lkeV, interpreted 
as a hot spot on the neutron star surface) becomes weaker until non-detection. In the simultaneous observation of the weak 
Chandra /RXTE spectrum, there is no more indication for the ~ 1 keV thermal component while we detect a colder thermal 
component (~0.4keV) that we interpret as the emission from the cold disc. A hint for a 6.4 keV iron line is detected, together 
with an excess around 6.8 keV and absorption feature around 7.1 keV. The latter two features have never been detected in the 
spectra of accretion-driven millisecond pulsars before and, if confirmed, would suggest the presence of an expanding hot corona 
with high outflow velocities. 



Key words. pulsars:individual IGR J00291+5934. 
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d 1 . Introduction tems supports the theory that LMXRBs are indeed the progen- 

itors of millisecond radio pulsars with a low magnetic field. 

In Low Mass X-Ray Binaries (LMXRB) hosting a neutron star, ^ AnN/rcn i a ^ ■ t 

J ' ^ Currently seven APMSPs are known. A r ecent review of 

the accreting star is thought to be very old and weakly magne- . u f a- vTr^ a »u 

° ° J b these systems can be found in Wimands (2005) . Among them, 

tised, compared to neutron star High Mass X-Ray Binaries. It ^.^^ .^^ ^ ^^^^ ^^.^ p^^.^^^ J0029 1+5934 is the fastest 

is believed that in LMXRBs, the accretion disc is generally not ^^^^^^ ^p^SP Evidence for the existence of rapidly (msec) 

influenced by the magnetic field and can extend very close to ^p^^^^^. ^^^^^^^ ^^^^ ^ ^.^^^ ^^^pj^ LMXRBs (at least 

the neutron star surface in a slow, "spot-less" accretion. This • ■ , , , . \ u u j j u »u j » 

^ in an additional 1 1 systems) has been provided by the detec- 
seems to be the case for the majority of LMXRBs for which ^ i u » n /"u * n "\ a ■ 

■' 1 tion of nearly coher ent oscillations ( bu rst osciUations ) during 

no regular pulsatio ns have been observed (see e.g. White et alj j^p^.j ^.^ay bursts ('chakrabartv"2004l). Nevertheless, only for 
[19951 | Psaltis| | 2004 for a review). In a few cases though, regu- ^ lMXRBs can we actually see the spin due to regular X- 
lar X-ray pulsations have been detected with spin periods rang- pulsations. The reason for this discrepancy is still debated 

ing from about 120 sec for GX 1+4 to less than 10 msec, e.g. ^ ■. ■ , , u u • » »u • • ^ *u v 

° ° and It IS not clear why the physics at the origin of the X-ray 

for SAX J1808.4-3658. The discovery of the latter accretion- , .. i i »u u»u r i 

■' pulsations reveals itself only through the presence of pulsations 

powered millisecond pulsar (APMSP) has been very important j » • »u * i j ■ * *u * 

^ F V / F and not in the spectral and tiimng properties of the sources that 

in the LMXRB evolution scenario: the existence of such sys- 



Send offprint requests to: A. Paizis: Ada.Paizis@obs.unige.ch 



' With the exception of the very recently discovered 
HETE JI900.1-2455 liMorean et al..2005iV 
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remain similar between APMSPs and non pulsating LMXRBs 
(Wiinands 2005). 

In this paper we report on observations of the fastest 
discovered APMSP, IGR J0 029 1+5934, performed with the 
Chandra X-ray Observatory dWeisskopf et alJl2002l) . 12 days 
after the discovery. To have a better overview of the spectral 
behaviour of the source prior to our observation, we have also 
analysed the available Ros si X-ray Timing Explorer (RXTE ; 
Llahoda & PCA Tea m 1996) observations, part of which were 
simultaneous to our ISksec Chandra observation. 

The paper is organised as follows: in section two we give an 
overview of the current knowledge of IGR J0029 1+5934 from 
its discovery. In section three we present our observations and 
data reduction methods. In section four we present our results 
that are then discussed in the last section. 

2. IGR J002gi +5934 

IGR J0029 1+5934 was discovered by the INTErnational 
Gamma-Ray Astr ophysics Laboratory (INTEGRAL; 
IWinkler et alJ l2003l) on December 2nd 2004, during rou- 
tine monitoring of the Galactic plane (JSckert et al.. .2004.) . 
Follow-up observations with RXTE r evealed the presence of 
coherent pulsations at ~598.89 Hz ('M arkwardt et alJ l2004 
with an energy dependent fractional amplitude. Further analy- 
sis and observations with RXTE showed that the neutron star 
is in a 8844 s (2.45 h) orbit. The upper limit of 0.16 Mq on 
the mass of the companion star in IGR J0029 1+5934 implies 
that the companion is most probably a hot brown dwarf 
jGallowa v et al.'2005V 

The outburst X-ray spectrum could be fitted with an ab- 
sorbed power-law with F = 1.8 a nd fixed column den- 
sity A?H=2xl02'cm-2 for INTEGRAL dShaw et al.l l2005') and 
F = 1.7 and measured c olumn density A^h=1-7x10^' cm"^ 
for RXTE dGallowav et al.ll2 005). More physical models (ther- 
mal Comptonisation) led to poorly constrained el ectron tem- 
peratures with a cut-off between 80-120 keV jShaw et al.l 
l2005l) . Llonker et alJ (l2005l) observed IGR J00291+5934 with 
Chandra one month after the discovery of the source, likely 
at its quiescent flux level, obtaining a 0.5-lOkeV flux of 
~10"^^ ergcm"^ s 'and a neutron star effective temperature of 
-0.3 keV. 

Inspection of the RXTE/ASM archive showed that the 
source was likely also active in the past on two occas ions, lead- 
ing to a tentative recurrence time of about 3 years jRemillardl 
12004 ). 

Observations with ground based telescopes allowed the 
counterparts at optical and radio wavelengths to be dis- 
cover ed (.Poolev .2004 iFender et al.. .2004.: Fox & Kulkarnil 
I2OO4I) . The most accurate optical position of the source is 
RA^OO** 29"* 03.0 6^ and Dec=+59° 34' 19.0" (0.5" uncer- 
tainty) (iFox & Kulkarni.2004) . 

3. Observations and data reduction 

An X-ray intensity history of the 2004 outburst of 
IGR J0029 1+5934 is presented in Fig. [T] where the 
INTEGRAL detection followed by the RXTE follow-up ai-e 
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Fig.l. X-ray intensity of IGR J0029 1+5934 throughout the 
2004 outburst. The open triangle indicates the discovery with 
INTEGRAL . The 2.5-25 keV PCA flux evolution is shown 
(open squares).The time of our Chandra observation is indi- 
cated by the arrow in the right. After MJD 53355 the source 
flux is strongly influenced by residual flux and variability of (at 
least) flie nearby V709 Gas. 

shown. The time of our Chandra observation is indicated by 
the arrow on the right. As can be seen, we observed the source 
at a very faint flux level, after the source had decayed signifi- 
cantly towards quiescence. 

3.1. The Chandra data 

The 18ksec Chandra observation of IGR J0029 1+5934 was 
performed on 2004 December 14th from 02:43:41 until 
08:02:38 UT (MJD 53353). We reduced the data in a stan- 
dard manner using the CI AO version 3.2 software package 
and CALDB version 3.0.0. The obtained spectra ha ve been 
analysed with the ISIS analysis system, version 1.2.4 llHouckl 
.2002) . 

We used the Hi gh Energy Transmis sion Grating 
Specti-ometer, HETGS fCan izares eTall 12000"). It has two 
sets of gratings, the High Energy Grating, HEG, and Medium 
Energy Grating, MEG, covering the energy ranges of 0.8- 
lOkeV and 0.4-5. OkeV, respectively. The focal plane imager 
used is the Advanced CCD Imaging Spectrometer (ACIS-S), 
an array of six CCD detectors normally used as readout for the 
photons dispersed by the gratings. The CCDs were operated in 
a sub-array mode where only half the CCD was read out. This 
did not affect the dispersed spectra, but served to reduce the 
frame integration time from the usual 3.2 s to 1.7 s, and thus 
reduced the presence of pileup (see below). 

Given the low signal to noise obtained, we extracted the 
zeroth-order (undispersed) spectrum and the first order dis- 
persed spectra {m = + 1 for HEG and MEG) for a total of five 
spectra. Higher order spectra were not considered. The zeroth- 
order image of the source is well resolved, but is also mildly 
piled up (pile-up fraction of 18%). In all spectral fits wherein 
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we fit the zeroth-order data, we account for these effects by 
using the ISIS implementation of the pileup model of Davis] 

Within the field of view of ACIS-S we detect an- 
other source, 17' away from the zeroth-order image of 
IGR J00291H-5934 . The source is V709 Cas, a cataclysmic 
variable. The zeroth-order image of V709 Cas lies on top of 
the HEG m = -1 dispersed spectrum of IGR J0029 1+5934. 
The position of V709 Cas corresponds to the l-1.3keV region 
of the HEG m = -1 spectrum, that has been ignored in the 
spectral analysis. Outside of 1-1.3 keV, the removal of events 
associated with V709 Cas is extremely efficient and no contam- 
ination is occurring. 

Inspection of the 5 spectra showed that the separate grating 
arms agree with each other This allowed us to merge the two 
HEG (m - +1) and MEG (m = +1) spectra into two combined 
spectra in order to increase the signal to noise ratio. This led to 
three final spectra that we used in the analysis (combined first 
order HEG, combined first order MEG and zeroth-order). For 
the fit, we binned the data to obtain a minimum of 16 counts 
per bin for both HEG and MEG as well as a minimum number 
of 16 channels per bin for HEG and 8 for MEG. 

The 1 .67 msec pulsations cannot be detected in our data, 
due to the 1.7 s frame integration time in our observations. We 
have used both 'unbinned techniques' for searching for vari- 
ability on all time scales > 1.7 s (e.g., the Bayesian Blocks 
method of Scars in prep.^) as well as 'binned' period 

folding techniques to search for variability on time scales com- 
parable to the orbital period. No statistically significant vari- 
abiUty was found. 

3.2. The RXTE data 

One day after the discovery of IGR J0029 1-1-5934, a regular 
monitoring of the source was performed with RXTE via Target 
of Opportunity (ToO) observations starting from December 3rd 
until December 16th. The RXTE observation on December 14th 
was performed simultaneously to our Chandra pointing (from 
03:08:48 to 08:53:52 UT). The RXTE data were extracted with 
the HE ASOFT software, ve rsion 5.3.1 using our standard pro- 
cedures JWilms et alJl999l) . 

The simultaneous Chandia/RXTE fit was performed in 
the already mentioned ISIS package. We rebinned further the 
HETGS spectra to obtain a minimum number of channels per 
bin equal to 64 for HEG and 32 for MEG. In this way the num- 
ber of PCA and HETGS data bins is comparable and during 
the fit similar weight is given to all the available spectra (HEG, 
MEG, PCA). We fitted the PCA spectra in the 3-13keV spec- 
tral range with a 0.5% systematics and grouped the data to have 
a minimum signal to noise of 5. The interpretation of the PCA 
December 14 spectrum of IGR J0029 1-1-5934 is dehcate due to 
the presence of (at least) the nearby cataclysmic variable V709 
Cas (17' away, see Section 3.1) that is clearly detected in the 
Chandra field of view. 



^ See the implementation of this algorithm in ISIS/s-Iang at 
^http://space.mit.edu/CXC/analysis/SITAR/_ 



To take into account the contamination from nearby 
sources, we subtracted the last RXTE pointing avail- 
able from all the previous observations. The integrated 
2-8 keV flux in this last PCA observation is consistent 
with the Chandra determined flux of V709 Cas alone 
(2.4x10"" erg cm"^ s"'), independently confirming the non- 
detection of IGR J0029 1-1-5934 in this observation. 

4. Results 

We extracted the X-ray position of IGR J00291H-5934 from 
the zeroth-order image, obtaining RA^OO*" 29™ 03.08' and 
Dec=H-59° 34' 19.2" equinox J2000 (90% confidence error 
of 0.6"^'). This position, compatible with the optical and ra- 
dio ones, was imm ediately announced to the community by 
iNowaket al.l(l2004 . 

The first order HEG and MEG spectra of 
IGR J0029 1+5934 are shown in Fig. H The best fit 
obtained, also shown in Fig. |2l is an absorbed power- 
law with column density"* A^h=(4.3+0.4)x10^' cm"^ and 
r = 2.06 + 0.07 with a reduced xl^0.19 for 431 d.o.f. The 
absorbed flux is 2.5x10"'' erg cm"^ s~' in 0.5-8keV and 
1.9x10"" erg cm-2 s"' in 2-8 keV (-ImCrab) with about 10% 
uncertainty. 




0.5 1 2 5 



Energy (keV) 

Fig. 2. Chandra HEG (order m = +1 combined, black) 
and MEG (order m - ±1 combined, blue) spectra of 
IGR J0029 1+5934 . The model shown is an absorbed power- 
law. 

Residuals in the continuum are found by taking the best- 
fit continuum model, and then using that as a 'Bayesian prior' 
for the expected count rate in each bin of the unbinned spec- 
trum. One can then apply the Bayesian Blocks algorithm of 
Scargle (2005, in prep.) to search for the most significant data 
residuals^. Using this procedure, we find the most significant 

^ See http://asc.harvard.edu/cal/ASPECT/celmon/index.html 
The improved model for the absorption of X-rays in the interstel- 
lar medium "tbvabs" Version 1.0 was used (Wilms et al. 2000) 

' This procedure is described in detail at: 
http://space.mit.edu/CXC/analysis/SITAR/bb-experiment.html, 
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residuals to be excess emission between approximately 6.4 and 
6.8 keV (1.94-1.82 A), and an emission deficit between 7.0 and 
7.2 keV (1.73-1.78 A), as shown in Fig.g] Fitting just the HEG 
spectrum between 1 .5 and 3 A, where we have grouped the data 
uniformly by 7 bins, and using the statistic of Cash ( 1979) , if 
we set the normalisation of the absorption line to zero and refit, 
the statistic increases by 10.8. This is approximately the 99.3% 
significance threshold. (Setting the normalisation of the emis- 
sion line to zero and refitting only changes the Cash statistic 
value by 3.) Thus, barring systematic uncertainties of the de- 
tector, this absorption feature is approximately 3cr significant. 
Hints for the presence of these featues are seen in the zeroth- 
order spectrum as well, as we show below. 




1.7 1.8 1.9 2 2.1 2.2 



Wavelength (A) 

Fig. 3. HEG data uniformly rebinned by a factor of 7, where 
we have fitted this region with a power law plus broad line plus 
absorption feature. The latter is the most significant feature in 
our HEG spectra (approximately 99.3% significant) and is also 
present in the zeroth-order spectrum. 

The simultaneous HETGS-PCA spectrum gave an overall 
harder spectrum with F = 1.8. The harder continuum set by 
PCA gave a soft excess in the Chandra spectrum, which was 
then adjusted in the fit resulting in a smaller column density 
(A^H=3.1 X 10^' cm"^). Nevertheless the fit was not very good, 
reduced of 1.35 (163 d.o.f.), and indeed the low end of the 
MEG spectrum is not well described; the new A^h is lower than 
the value obtained fitting the Chandra data alone and predicts 
more photons in the soft end of the spectrum than what actually 
detected by the MEG. Furthermore, the higher end of the PCA 
spectrum shows a slight excess above 8 keV, meaning that the 
1.8 slope is not enough for the PCA data. 

To adjust for these issues, we added a soft the rmal part to 
the spectrum (DISKBB model. lMitsuda et alJl984 obtaining a 
very good description of the data (j(;l-0.9, 161 d.o.f.): the PCA 
data are better described by the new power-law (r= 1.7) and the 
softer part of the spectrum is well fitted by A'h=4.3 x 10^' cm"^ 
and a thermal component of kT - 0.42 keV (13% contribution 
to the absorbed 0.5-8 keV flux) that compensates for the soft 



excess induced by the hard slope. Fig.l^shows the HEG, MEG 
and PCA specti-a of IGR J00291-H5934 together with the best 
fit model. 

The discrete feature at 2. 1 keV is a known calibration origin 
and is due to a jump in the HETGS effective area. The excess 
at ~5 keV has no known instrumental origin and could be real; 
however, it is seen neither in the HEG nor in the zeroth-order 
spectrum and is less significant than the akeady discussed iron 
region features that are mirrored also in the zeroth-order spec- 
trum. 




0.5 1 2 5 10 



Energy (keV) 

Fig. 4. HEG (black), MEG (blue) and PCA (red) spectra of 
IGR J0029 1-1-5934 . The fit shown is a composition of an ab- 
sorbed power-law and disc blackbody. 



We find that the normalisation of the PCA spectrum is 
about 13% higher than for Chandra . Between these two mis- 
sions, normalisation factor up to 20% have been found before 
(e.g., Juett et al. 2003, and references therein). We note here 
that if we use a non-cleaned PCA spectrum, i.e. not corrected 
for the contribution of the nearby sources, we obtain a normal- 
isation factor of 50% with respect to Chandra. This high value 
is most likely due to the contamination from nearby sources 
and background in the PCA spectrum. 

As a last step, we included the Chandra zeroth-order spec- 
trum that is affected by pile-up and thus the spectral shape from 
the source is distorted. Fig.|5lshows the HEG, MEG, PCA and 
zeroth-order count spectra. The spectral distortion due to the 
pile-up in the zeroth-order spectrum is visible (dip in the 1- 
2keV region in the green spectrum). The emission/absorption 
features detected by the HEG in the iron region (already dis- 
cussed. Fig. O are mirrored in the zeroth-order as well. The 
information that is lost due to pile-up can be partly recov- 
ered using kernel models that "correct" the distortions IPavisI 
C200ll). Nevertheless, the result cannot be perfect and the in- 
clusion of the zeroth-order spectrum can make the overall fit 
quality slightly worse. We obtain an absorbed power-law and 
thermal component, compatible with the previous case, with 
column density A?H=3.9x 10^' cm^^^ p ^ j 7^ ^7^.^ ^ 0.4keV 
and a reduced ;i^y= 1.1 5 for 293 d.o.f. Again we find a PCA 
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normalisation of about 13% higher than for Chandra. Table ^ 
summarises the results of our spectral analysis. 




Energy (keV) 

Fig. 5. HEG (black), MEG (blue), PCA (red) and HETGS 
zeroth-order (green) count spectra of IGR J0029 1+5934 . Note 
the emission/absorption line features already discussed (Fig.|3 
in both the zeroth-order and HEG spectra. The break in the 
HEG spectrum (data points between 1-1.3 keV) is due to the 
fact that in this range we have used HEG m - +\ alone and not 
the averaged m - +1 orders, due to the contamination by V709 
Cas in the HEG m - -1 order. 



Table 1. Best fit parameters for the absorbed power-law model 
of the HETGS spectra of IGR J00291-H5934 and for the ab- 
sorbed power-law plus thermal component of the HETGS, 
PCA and Chandra zeroth-order spectra. A^h is the column den- 
sity in units of 10^^ cm~^; ^r,„ is the inner temperature of the 
thermal (DISKBB) component; d.o.f. = degrees of freedom 
(the decrease in d. o. f. in the combined HEG / MEG and PCA 
spectrum is due to the heavy rebinning we performed in the 
HETGS data, see text); in all the fits we fixed to 1 HEG, MEG 
and zeroth-order normalisation factors. The indicated errors are 
at Icr. 



HEG / MEG 



HEG / MEG / PCA 



A'h 4.3±0.4 

r 2.06±0.07 

PCA norm 

reduced xl 0.79 (431 d.o.f.) 



4 3+u./ 

^•-^-0.5 

1.7±0.08 

0.42!°;°', keV 
1.13 

0.9(161 d.o.f) 



HEG / MEG 
PCA / 0th order 



A'h 

r 

kT,„ 

PCA norm 
reduced 



n q +U.3 

■'•^-0.4 
1.7±0.04 

0.40±0.06keV 

1.13 

1.15 (293 d.o.f.) 



To have an overview of the outburst behaviour of 
IGR J0029 1+5934 before our Chandra observation, we anal- 
ysed all the available RXTE/PCA observations. The PCA spec- 
tra of IGR J0029 1+5934 can be fitted by a combination of a 
thermal component (~1 keV) and power-law (r~1.6)^. A hint 
for an iron line at 6.4 keV (fixed) is detected. During the out- 
burst, while the source is decaying, the power-law photon in- 
dex shows no particular trend (unlike the normalisation that 
decreases significantly) while the thermal component becomes 
weaker and can be constrained no more by PCA. Already two 
days before our Chandra observation the RXTE observations 
are consistent with a single power-law with a marginal detec- 
tion of the 6.4keV line. We believe that in our Chandra spec- 
trum we still see a residual of the disappearing 6.4 keV line as 
well as a soft component (0.4 keV) that cannot be constrained 
by PCA alone. 



5. Discussion 

We have studied the spectral evolution of the accretion- 
powered millisecond pulsar IGR J0029 1+5934 with 
Chandra and RXTE during the December 2004 outburst 
that led to its discovery by INTEGRAL. 



* This slope is adjusted to the softer 1.7 value when fitted with the 
Chandra simultaneous spectrum of December 14. 



5.1. Long term variability 

At its discovery IGR J002 9 1+5934 had a 5-lOOkeV flux 
of about 10" erg cm -2 s-' (,ShawetalJl2005l) that decayed to 
about 10 "'erg cm"^ s '(extrapolated from our HETGS/PCA 
best fit) about 12 day s later. Assuming a distance of 5kpc 
dCiallowav et all l2005l) we obtain that the X-ray luminosity 
changed from Lx~10^*ergs 'to Lx~10^^ erg s '. This corre- 
sponds to a change from ~0.01 to ~0.001 LEdd that places 
IGR J0029 1+5934 at the lower end of the dim AtoU sources 
(-0.01-0.3 LEdd), similarly to other APMSPs. Such a peak lu- 
minosity is very low if we consider the bright neutron star 
X-ray transients (XRTs) that display outbursts with peak lu- 
minosities of Lx~10^'-10^^ ergs"'. The history lightcurve of 
IGR J0029 1+5934 from the RXTE/ ASM indicates that the flux 
measured at the time of the discovery is indeed the peak of the 
outburst but unfortunately, the distance to IGR J0029 1+5934 is 
not well constrained and we cannot be sure of the absolute 
value. It may seem reasonable to expect that the peak lumi- 
nosity is indeed low: APMSPs are very compact systems, the 
disc outer radius R is relatively small compared to wider or- 
bit LMXRBs and this limits the total disc mass (ocR^) that can 
build up in quiescence. When the outbur st is triggered, the mass 
accretion onto the central object is ocR^ ( lGierhnskietal.l2002l 
and references therein). Nevertheless, things are not so straight- 
forward and there are cases where the above scenario does not 
hold, e.g. in the case of the ultra compact (1 1.4 minute period) 
LMXRB 4U 1820-30 that can still reach X-ray luminosities 
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highe r than a few times 10^^ erg s ' jBallantvne & StrohmaveJ 
l2004 . 

APMSPs differ from the bright XRTs also in the 
quiescent luminosity: standard XRTs range between 
Lx-lO^^-lO^^ergs^' (Campana et al. 2002) while APMSPs 
seem to be dimmer than ~10^^ erg s^' dWiinan ds 2005). Our 
Chandra measurement of IGR J0029 1 +5934 gives a 0.5-8 keV 
flux of ~10"" erg cm"^ s"' that, com pared to the 0. 5-lOkeV 
flux of ~10 ergcm"^ s 'by Jon ker et alJ ( l2005l) . clearly 
shows that we detected IGR J0029 1+5934 during its outburst 
decay, prior to its quiescent state. 

5.2. The spectral continuum 

The study of the other APMSPs has shown that during the 
outburst the spectra can be well fitted by a two component 
model, a soft thermal component and a harder one. In gen- 
eral the soft component is interpreted as blackbody emis- 
sion from a heated hot spot on the neutron star, responsi- 
ble for the X-ray pulsations, and ranges within about 0.6- 
1 keV (see e.g. G ierlihski et al. 2002; Gierliriski & Poutanen 
l2005UTuett et'an i2003). The harder component is interpreted 
as thermal Comptonization by a plasma heated by the accre- 
tion shock as the material, collimated by the magnetic field, 
impacts onto the neutron star surface. The source of seed pho- 
tons for Comptonisation seems to be either the hot spot it- 
self ( Gierlihski & Poutanen 2005) or the colder accretion disc 
Jritarchuk et al.l2002ft . 

In our analysis, we find that the spectra of 
IGR J0029 1+5934 can be fitted by a combination of a 
thermal component and a power-law. Along the outburst, while 
the source is decaying, the power-law photon index shows no 
particular trend while the thermal component becomes weaker 
and can be constrained no more by PCA. At this point, we 
obtained our Chandra observation and the simultaneous weak 
Chandia/RXTE spectrum can be well described by the combi- 
nation of a colder thermal component (kTi,, = 0.4 ke V instead 
of ~1 keV) and a power-law. The presence of the power-law 
means that if there is a cut-off in the spectrum, then this occurs 
at higher energies than the HETGS/PCA range and we are not 
able to constrain it. A Comptonisation model in our data gives a 
poorly constrained Comptonising plasma temperature of about 
50keV. The presence of the thermal component means that 
not all the available soft photons go through the Comptonising 
medium, instead there are parts that are seen directly, sirnilarly 
to wha t w as found in other APMSPs in outburst, iMiller et alJ 
("2003") and'Gierlihski & PoutanenI (l2005l) on XTE J 1 75 1 -305, 
[Gierhhski et al. (2002) on SAX J1808.4-3658 and Juett et al. 
I2OO3I) on XTE J0929-3 14. 

We interpret the evolution we see in PCA data and the final 
simultaneous PCA/Chandra spectral properties in the follow- 
ing way: at the on-set of the outburst the accreting matter is 
channelled onto the neutron star surface and a hot spot is cre- 
ated. We have indications for the hot spot in the early PCA data 
with a thermal component around 1 ke V that is likelv to be re - 
sponsible for the pulsations detected bv lGallowav et al.^ (20o3). 
The disc is colder (<1 keV) and we do not detect it dii^ectly in 



the PCA spectra alone, but we see its effect in a hint for the 
6.4 keV line in the PCA data, most likely originating from irra- 
diation of the cold accretion disc by the X-ray source.^ Along 
the outburst, the hot spot becomes progressively weaker and al- 
ready 10 days after the on-set of the outburst we are not able 
to const rain it with PCA anym ore. This is consistent with the 
fact that lGallowav et all (l2005i) detect no pulsations in the last 
phases of the outburst. The disc is likely becoming colder as 
it is illuminated by the fading X-ray source and is possibly re- 
ceding after the outburst, similarly to other LMXRBs in which 
at a low accretion rate, the disc is more distant from the com- 
pact object and the overall spectrum is hard. At the time of our 
Chandra observation the hot spot (along with its pulsations) 
is most likely off and we detect a cold accretion disc around 

0. 4 keV. 

5.3. Discrete features in the spectrum 

Despite the regular search for discrete features in APMSPs 
in the past, no line was found besides the 6.4 keV one in 
SAX J 1808.4-3658 jGierlihski et alJ2002h . In our study, a hint 
for an iron Ka line (fixed at 6.4 keV) is detected in the PCA 
observations, while a marginally significant excess in the 6.4- 
6.8 keV region and a 3(T significance deficit around 7.1 keV are 
detected in both the HEG and zeroth-order Chandra spectra 
(Figs.|3]and|5j. Emission lines at 6.4 keV have been observed 
from many X-ray pulsars and are believed to be produced by 
fluorescence of weakly ionised iron (less than Fe XV). A likely 
candidate for this cold matter is the accretion disc that can be 
also responsible for the iron K-edge ~7.1keV feature we de- 
tect in the spectrum. The emission line that we detect around 
6.8 keV (if real) would require the presence of a highly ionised, 

1. e. hot, corona in IGR J0029 1+5934, consistent with the high 
plasma temperature (>10keV) inferred from the spectral study 
performed here, as well as in Shaw et al. (2005). This scenario 
would require a different medium for each feature of Fig. 
(cold disc for the 7. 1 ke V feature and hot corona for the 6.8 keV 
one). 

Another possible, although speculative, interpretation of 
the discrete features in Fig.|3lcould be the following: the emis- 
sion feature we detect around 6.8 keV could be the redshifted 
part of the 6.97 keV line (expected from Fe XXVI) while the 
absorption feature around 7.1keV could be the blueshifted 
6.97 keV line. In this scenario we would basically obtain that 
the 6.4 keV line is produced by the cold disc while the 6.8 
and 7.1 keV lines are a P-Cygni profile from an expanding hot 
corona with outflow velocities of about 6000 km/sec. Such (and 
higher) outflow velocities are possible and have already been 
observed in cataclysmic variables (e.g. .'Woods et alJl992h and 
in the X-ray binary SS433 (Migliari et al. 2005). The P-Cygni 
profile scenario is consistent with all the available data we cur- 
rently have on IGR J0029 1+5934. However, the relatively low 

^ A direct simultaneous measure of the hot spot and disc has 
been reported in the case of the APMSP XTE J1751-305 by 
Icierlinski & Poutanen ( 200^ where two soft components were de- 
tected in the broad-band XMM-Newton/RXTE spectrum: -0.6 keV, as- 
sociated to the disc and ~ 1 keV, the hot spot. 
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statistical significance of the discrete features does not allow us 
to firmly establish this interpretation and underUnes the impor- 
tance of obtaining quick follow-ups of these transient events. 
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